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Abstract. The Noether gauge symmetry approach is revisited to study various 
quintom scenarios (those that arise by the presence of two dynamical scalar helds) 
to comprehend the role of dark energy in our universe. For such models, we obtain 
smooth parameterizations of the equation of state of dark energy across the boundary 
of cosmological constant wa = —1- This study gives rise to two new cases of the 
potential V (<(), a ), due to a quintom field in which nonlinear coupling of the scalar 
fields arise. Besides we report that a few cases of Noether gauge symmetries and their 
invariants in [Adnan Aslam, et. ah, Astrophys Space Sci (2013), 348:533-540] are 
incorrect. Consequently, the given cosmological model in their paper is not a feasible 
quintom model. 


1. Introduction 

Noether symmetries play an essential role in finding conservation laws (or invariants) of 
any physical system nsing Noether’s first theorem mm- For example, invariance of a 
Lagrangian of a physical system (in which force is dne to the potential in the field, i.e., 
conservative systems) nnder time translation and rotation correspond to conservation 
of energy and angnlar momentnm in favorable bonndary conditions. Noether symmetry 
approach provides an elegant way to find conserved qnantities for a given Lagrangian. 
Noether symmetries for a given Lagrangian generate possible forms of the nnknown 
fnnction of the field variables which help in redncing the dynamics and to find exact 
solntions. In [3], the anthors have nsed Noether symmetry approach to get possible 
cosmological solntions in /(i?)—gravity. On the other hand, Noether gange symmetry 
approach has an edge over standard Noether symmetry approach as it nnveils more 
conserved qnantities than from the other approach. This is clear from the observation 
that certain conserved qnantities do not correspond to invariance of Lagrangian nnder 
continnous gronps of transformations bnt nnder gange transformations. For example, 
conservation of charge in qnantnm field theory correspond to a gange invariance of the 
given Lagrangian. 
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An important problem in modern cosmology is to comprehend the role of dark 
energy that is now snpported by observational evidence that at present the universe is 
going through cosmic acceleration at a fast pace. There are three important evidences 
to support above argument based on the experimental study of (a) supernovae la 
(SNIa) [HI Ho] (b) cosmic microwave background radiation along with large-scale 
structure surveys (CMB & LSS) |8] and hnally (c) the age of our universe calculated by 
incorporating dark energy [2|. Einstein’s general relativity at large cosmological scales 
requires that our universe must carry two special features, i.e., it is both homogeneous 
and isotropic. Therefore, our universe has a mysterious form of matter whose role is 
opposite of gravity such that its pressure is negative that accounts for this acceleration. 
Observations conhrm the energy density of dark energy occupies 70% of our universe. 
It is also clear that the role of this mysterious matter at an earlier cosmological epoch is 
opposite to the present epochs because otherwise it may not have been possible to see 
present large-scale structures of our universe. 

Equation of state wde = Pde/pde, where pde and pde is the pressure and 
mass density, characterizes dark energy. For example, the case of a nonzero and 
positive cosmological constant boundary corresponds to wa = —1, in which case p\ 
is independent of the scale factor a{t). A quintessence held is a dynamical held for 
which evolution of its equation of state evolves in the range wq > —1. Similarly for a 
phantom held, wp < —1. There are evidences that the cosmological boundary is crossed 
[7] therefore one would be interested in knowing a dynamical held for which equation 
of state is able to evolve the cosmological constant boundary wa = — 1- On the other 
hand, there exists a no-go theorem that forbids the equation of state of a single scalar 
held to cross over the cosmological constant boundary [13] . One possible solution to this 
problem is to introduce a combination of two dynamical scalar helds, i.e., a canonical 
held (j) and a phantom held a. Such phenomenological models are known as quintom 
models which give rise to quintom cosmology in which the dynamics of dark energy is 
understood by using parameterizations of its equation of state [aisiE]. 

In this paper, our purpose is twofold. First we point out all the errors in [1] in the 
next section. Secondly we discuss some more cases in the third section that arise from 
Noether gauge symmetry approach in which the dynamical evolution of combination 
of the two helds (in the form of nonlinear sum or diherence) yield interesting quintom 
models. Later we discuss useful insights in quintom cosmology obtained from those new 
solutions. 

2. Basic Quintom Model 

For a simple double held (0, a) quintom model with an arbitrary potential 17(0, cr) [HIS], 
the action integral has the form 

S' = ^ y [R + 0;^0’'" - - 2V (0, cr)], 


( 1 ) 
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where g denotes the determinant of the Riemannian metric of spatially flat FRW metric 
g^^i, = diag (1, —a(t)^ ^ 3 ), and where metric of the three dimensional Euclidean space 
is ^ 3 . Further, a(t) is the usual scale factor describing the expansion of our universe. 
The field variables 0(f) and a{t) correspond to dynamical canonical and phantom fields, 
respectively [2]. It is important to include all physical equations in our analysis so the 
equations of motions are given by 

d /a V 

2 - + - = -p, 


( 2 ) 

(3) 

(4) 

(5) 


• dV 

0 + ?>H(f) H——— = 0 , 
d 0 

dR 

d + dHa - — = 0 . 

dcr 

The effective energy density and effective pressure are 

1 -,0 1 


P = 


P= 2 '#-' 


/)^ - + R( 0 ,(t), 

1 


2 ^ 


R(0,a). 

Thus we obtain the equation of state 


— a 


w = 


02 - d 2 + 2 R 

The necessary condition for any feasible quintom model requires that close to the 
cosmological constant boundary w = — 1 , we must have w'\w=-i 7 ^ 0 , in which case the 
dynamical evolution of both helds generates a quintom scenario such that the boundary 
is crossed smoothly. Due to its signihcance we obtain a general form of w', using above 
dehnitions of effective pressure P and effective energy density p, along with equations 
of motions ([2]). By differentiating ([5]), with respect to t, we get 

8(00 - dd)R - 4(02 _ + 0 R^) 


w = 


(02 - d2 + 2R)2 

which upon using equations of motions (| 2 ]), becomes 

-24RFr(02 _ ^2) _ 4(^2 _ ^2 ^ 2R)(010i + dl4) 


w = 


( 6 ) 


(7) 


(02 - d2 + 2R)2 

The above equation can be used to keep a consistency check on the cases especially 
when closed form of the solutions of equations of motions can not be obtained. We now 
discuss all possible cases in which equations of motions can be integrated explicitly with 
the help of Noether gauge conditions. 

In [1], the authors have used the Noether gauge symmetry approach to hnd Noether 
symmetries of the Lagrangian L(a, d, 0, 0, a, d) of ([T]), given by 

'1 -.o 1 . 


L = 


-3(2(2 T cP‘ 




v"( 0 ,ct; 


( 8 ) 
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The main objective of the paper was to obtain possible forms of the potential V (0, a), of 
a quintom held which are allowed by the presence of Noether symmetries. The Noether 
gauge symmetries are dehned in terms of the vector helds 


d ^ d d 

^ + “a; +’ 

provided they satisfy the Noether gauge symmetry conditions 


(9) 


+ LiDtT) = DtG, 


( 10 ) 


where G{t,a,(j),a), is an unknown function called the gauge function of the 
corresponding symmetry. Using Noether theorem the conserved quantity can be 
calculated from the formula 


I 


TL + 


x^T) 


dL 

dxf^ 


-G, 


( 11 ) 


where x^ = (a, 0, a) and = (T,a,(3,a) denote the variables of Lagrangian ([H]) and 
coefficients of the Noether symmetry ([9]), respectively. Note that Dtl = 0, upon using 
the equations of motions, where 


Dt = 


dt 


. d 

® o-^ 

oa 


d d 

’a^ + '^a? 


( 12 ) 


We hrst obtain the determining equations and give a list of the cases of Noether 
symmetries from [T]. We check the validity of each of their cases and show that in 
some cases the equations of motions are not satished. 


2.1. List of Cases 

The equation flTOl) yields a system of linear partial differential equations (PDEs) to be 
solved simultaneously. This is system (10) in [1], which we include here for completeness 


o' 

o"' 

o"' 

(13) 

0fT 'l(j) 0, 

(14) 

6a^ - = 0, 

(15) 

-6Q!0 + O^fla = 0, 

(16) 

o"' 

b 

+ 

CO 

( 17 ) 

of" (3t — G^ — 0, 

(18) 

Qac^t T Ga — 0, 

(19) 

ex T 2,aQ.a — a'Tt — 0, 

(20) 

3a + 2000 — aTt = 0, 

(21) 

3a + 2a7o- — aTt = 0, 

(22) 

3a^V a + 0^010, + + a?VTt + Gj = 0. 

(23) 


The above system was solved for different cases in [T]. We show that the cases (4) and 
(5) given in their paper, are incorrect while case (2) is incomplete. Further the given 
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potential (20) in Section 5 (Cosmography) [T] fails to satisfy the equations of motions. 
Below we list these cases in detail. 


1. In case (4), they have obtained a Noether symmetry 



(24) 


corresponding to the gauge function and potential 


G — Cit + C2, 

V{(j),a) = - Ci 0 + F(cr). 


(25) 

(26) 


From equation (j9]), we hnd that for this symmetry T = 1, a = 0, /? = 1/a^, 7 = 0, which 


clearly do not satisfy equation (fT6|) of the determining equations of Noether symmetries 


because = 0 but (3a = —3/a"^, which can not be zero and leads to a contradiction. 
2. Similarly the two Noether symmetries in case (5) are 



(27) 


(28) 


corresponding to the gauge function and potential 


(29) 

(30) 


G — Ci(j) + C 2 , 
!/((/), a) = F{a). 


If we choose G of the above form then equation (1T8|) . implies a^/3t — Ci = 0, which is not 
identically satished by the hrst symmetry Xi, for which (3 = 0, therefore Ci = 0 which 
makes the gauge independent of 0. Similarly, for X 2 , we have (3 = [cit/a? + F{a)) and 


a = 0, whereas equation ffT6|l implies that (3 can no longer be a function of parameter 
a. Therefore both cases are incorrect and the resulting operators are not Noether gauge 
symmetries of the corresponding Lagrangian. 

3. The case (2), is not complete as it contains two arbitrary constants whose 
ranges are not specihed. Furthermore possible forms of the arbitrary function are not 
characterized 



(31) 


For example, if we choose Ci = 0, C 2 = 1,C3 = ±1, and the function F takes the form 
given below then the Lagrangian has an extra Noether symmetry. 

4. In Section 4 (Exact solutions), authors have used case (4) with ci = 0, to obtain 
an invariant of the Noether symmetry. They have also given Euler-Lagrange equations. 
This choice of the value of constant makes the gauge function G = C 2 , and potential 
function as E = F{a). As we have proved that the given Noether symmetry for this 
case in [1], is not a Noether symmetry therefore the resulting invariant (equation (18) 
of their paper), namely 



( 32 ) 


Noether symmetry in quintom cosmology 


6 


is not an invariant because Dtl, on using the Euler-Lagrange equations 
4aa + + 2a^F{a) — 2a^ = 0, 

= 0 , 

dF{a) 


3a • 


a 

3a 


a — 


-a 


a da 

is a non-zero factor given by 


= 0 , 


(33) 

(34) 

(35) 


D,I 


3d(j) 

a 


(36) 


This is zero provided 0 = constant. They have obtained two exact solutions for this 
value of 0, with further choices on F{a) = 2 and F{a) = 0, both of which become 
sub cases of case ( 6 ). The case F{a) = c ^ 0, gives four Noether symmetries whereas 
F{a) = 0, gives hve Noether symmetries 


d ^ d d 


'I- i9 (9 

9^3 = + (p-rr ’ 

ocp oa 


= 


(37) 

(38) 


A = A 

(90’ ^ da ' 

It is a case when the potential V{(j),a), is zero. For completeness we include their 
invariants 


7i = I (ba^ - a^((9^ - cr^)j , (39) 

h = ^ (Qtd^ — 4aa — ta^{^‘^ — d^) j , (40) 

l3 = a^{(j)a - a^), (41) 

J 4 = 0^(9, (42) 

I 5 = -a^a, (43) 


and discuss possible exact solutions as well as their cosmological implications in Section 
4. Similarly for F ^ 0, there are four Noether symmetries and 


correspondingly there are four invariants 

Ji = |(6a2-a2((92-d2 + 2C')) , (44) 

1 3 = a^{(j)a - a(j)), (45) 

1 4 = (46) 

h = -a'd. (47) 


5. In Section 5 (Cosmography), the only potential taken in their paper was expected 
to reveal some new results but it contains serious flaws. The potential was conhned to 
have the form 


1 /( 0 , a) 



0 )(a + 0 -l- 1 ), 


(48) 


with some specihc choice of constants. The above function has the derivatives 


d 0 




( 49 ) 
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therefore, in their paper, equations of motions ( 21 ) and ( 22 ) which come from equations 
(3) and (4) are wrong. Consequently, equations (24) and (25) are also incorrect. Hence 
the subsequent analysis is totally unnecessary and contains fallacies in all diagrams and 
cosmological implications in Section 6 (Stability). 

3. Additional Cases 

The system of determining equations flT^ — fl2^ . is difficult to solve in general. We have 
used Computer Algebra System (CAS) Maple—17, to carry out case-splitting with the 
help of an important algorithm ‘rifsimp’ that is essentially an extension of the Gaussian 
elimination and Groebner basis algorithms that is used to simplify overdetermined 
systems of polynomially nonlinear PDEs or ODEs and inequalities and bring them 
into a useful form. It turns out that the system is inconsistent with the choice of a time 
dependent gauge function. We have also verihed all the cases in the previous section. 
Furthermore we obtain some intriguing forms of the potential function with the use of 
GAS that can shed light on the physical insights of canonical and phantom helds and 
their behaviour at different epochs of our universe. To the author’s knowledge these 
results of dynamical invariants have not been reported before. Below we list these new 
cases. 

1. We now hnd a solution of system flT^ — fl2^ . in which a time varying invariant 
is obtained for a new quintom model. It is the case of a constant gauge function and 
potential V (</>, a) with a nonlinear power of sum of the canonical held ( 0 ) and phantom 
held (cr), i.e., (j) + a. Since we obtain a similar result in the case of diherence of both 
helds (0 — a) with a change of plus/minus sign, therefore we include both in a unihed 
form given below 


(50) 

(51) 


G = c, 

r('^,ff) = (i^±ff)”, riyto, 


where n is a non-zero integer. It would be interesting to see the impact of above potential 


due to the presence of both canonical and phantom helds on the density and pressure of 
the universe [2]. The resulting model is such that the sum (or diherence) of contribution 
of canonical and phantom helds on the potential is increasing for n > 0 , and decreasing 
for n < 0. Since the helds are time-varying in general therefore it could be that their 
sum (or diherence) is increasing for certain periods of time and decreasing for other 
periods of time. The value n = 0, corresponds to the case ( 6 ) in [1], that contains 
a typo, namely the Noether symmetry = F{a)d^, is incorrect as it contain F{a) 
which is not possible. It is X 4 = d^. We now obtain three Noether symmetries of the 
Lagrangian ([ 8 ]) for the above forms fl50|) & fl5T|) of the gauge and potential functions 


which are 



50 da 



( 52 ) 
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where these satisfy the commutator relations 

[Xi, X2] = 0 , [Xi, X3] = Xi, [X2, X3] = ±^X2. ( 53 ) 

The Noether symmetry Xi corresponds to translations in time whereas X 2 is the 
difference (and respectively sum) of the translations in the held variables 0 and a. 
Interestingly, X 3 is a non-trivial Noether symmetry which includes a mixture of both 
scaling (in t and a) and a hyperbolic rotation (in 0 and a). 

Noether theorem ensures that there exist three invariants (conserved quantities) 
corresponding to three Noether symmetries for the Lagrangian (j 8 ]) 


L 


—3aa^ + 





(54) 


From equation ffTTl) . the three invariants relative to each Xi,X 2 and X 3 of the above 
Lagrangian are 


h 


3aa? + 



(0±ct) 


(55) 



J 3 = 2 a^a — 3tad? + 




(56) 

(57) 


The hrst two invariants could be guessed from the form of the Lagrangian fl5T)) . However 
the third invariant is a non-trivial dynamical invariant and represents a quantity which 
changes over time but remains conserved during evolution. Therefore both Ji and I 2 
are hrst integrals whereas is a constant of motion. 

2. We also obtain another case in which the system has a solution with the gauge 
and potential functions of the form 


G = c, 

(6 zf aN 

^n^0,m^0, n ^ -m 

(0 ± a)™ 


(58) 

(59) 


where n and m are two positive integers, and we again write down two separate cases 
in a unihed form. The behavior of helds is more interesting in this case, for example 
if u = m, then V is the ratio of the difference to the sum of the helds (or in the other 
case it is the ratio of the sum to the diherence of the helds). For n ^ m, we obtain 
two distinct cases. If u > m, then the impact on potential of diherence of the helds is 
larger than their sum (or impact due to sum of the helds is larger than their diherence 
in the other case). On the other hand if n < m, then the impact of diherence of the 
helds is smaller than their sum (or impact due to sum of the helds is smaller than their 
diherence in the other case). The two Noether symmetries of the Lagrangian ([H]) for the 
above forms (l58|) & (l59|) of the gauge and potential functions are 




d d ad 2 f d < d \ 

df ^ ^ dt~^ 3 da n + m\ dcf)^ da J ’ 


( 60 ) 
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where the algebra is closed 

[Xi,X2] = Xi. (61) 

The two invariants corresponding to two Noether symmetries for the Lagrangian ([8]) 

(62) 


T Q-2| 3/^12 1-2 (0T<^) 

L = -3aa +a ( -0 - -a - —-—- 

2 (0 ± cr) 


have the form 


Ji = 3aa^ + 


— 0^ (0 =F cry 


(0±a) 


I 2 = 2a^a — 3taay + 


a? I t 


(0Ttr)^ 

(0±cr)' 


+ 




± 


n + m 


a — a 


(63) 


(64) 


In this case Ji is the hrst integral whereas I 2 is a constant of motion. 


4. Quintom Solutions 


Case 1. (1/(0, a) = 0) 

This is the case in the absence of a potential held, i.e., 1 /( 0 , a) = 0, in which we obtain 
hve invariants corresponding to hve Noether symmetries possessed by the Lagrangian 
([8]). We include this case to verify previous cases as well to comprehend the role of 
invariants Jj, and constant gauges that result the evolution of both helds. There are 
hve invariants for the case of constant gauge function, therefore we write Jj = q, where 
i = We use small letters to denote constancy of invariants and integration 

constants are represented by capital letters. Note that these constants may carry any 
constant value nevertheless non-trivial solutions of equations of motions depend if these 
constants are zero or not. For example, the invariant /4 = C4 , gives 0 = constant, 
if C4 = 0 and 0 = 04 / 0 ^ otherwise. Their values are specihed by certain physical 
constraints. To obtain all possible cases exhaustively we have carried out case splitting 
using computer algebra system and found following solutions, 
la. (i) (cj y 0 for all i) 

a{t) = ^ (6cit^ - 4c2t + , (65) 


a{t) = 


4ciC5 


\/Q{cl-cl) 


arctan 


3cit — C 2 

V6(ci-c|) 


+ C 2 , 


( 66 ) 


0(f) = 


4 C 4 C 1 


;arctan 


3cit — C 2 


^/Q{ci-cl) \VH4-yi) 
cl + 6(c2 - cl) 


C 3 — C 4 C 2 
C 5 


where Ci = 


12ci 


C-S ^ C 4 


In the absence of potential {V) the dynamical w reduces to a constant value 




w = — 


/)2 _ (j2 


= 1 


(67) 

( 68 ) 


(69) 
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Thus the simplest case yields p = p, which corresponds to a radiation model and the 
evolution of both helds is determined by the forms above. Subsequently in all sub-cases 
of case 1, the equation of state possesses the same form, 
la. (i) (cj 7^ 0 for all i) 


a(t) 

a{t) 


i(6cl^2-4c2^ + C'l)'/^ 

4ciC5 I 3cit - C2 \ 

- , = arctan — , = 


+ C 2 -, 


(70) 

(71) 


0(t) = 


4c4Ci 


v/6(ci-c|) 


arctan 


1 ry ^2 + 6(c| 

where Ci = — ^ 


3cit C2 I 

y6(ci-ci) 

—, Icsl > |C4 


C3 - CiC2 
C5 


12ci 

In the absence of potential (V) the dynamical w reduces to a constant value 


w = 


f )2 _ (j2 


= 1 . 


(72) 

(73) 

(74) 


Thus the simplest case yields p = p, which corresponds to a radiation model and the 
evolution of both helds is determined by the forms above. Subsequently in all sub-cases 
of case 1, the equation of state possesses the same form, 
la. (ii) (cj 7^ 0 for all i) 


a(t) 

a{t) 

4>{t) 


= i(6cii"-4c2i + 8Ci)‘'“, 


- 1 

' 3cit - C2 1 

1 + C2, 

Mci - 4) ' 

iy V6(cl - ci) ) 


^ 3cit - C2 \ 

1 , C3 - C4C2 

V6(cl - ci) ' 

\^Q{cl-cl)) 

C5 

1 ^7 *^2 ~ 6(c| - 

where Ci = ^ 

^ 1 0 _ 

1-^. |C4|>|C5| 


(75) 

(76) 

(77) 

(78) 


lb. (cj 7 ^ 0 for all i = 2,.., 5, Ci = 0) 

«(^) = ^ (Cl - 4c2f)^C ^ (79) 

(T(t) =-—In (Cl - C 2 t)-h C 2 , (80) 

C 2 

(f){t) = C2t)+ — —(81) 

C 2 C 5 

where C 2 = ±\/cl - c|, |c 4 | > \c^\ (82) 


which yields two solutions relative to positive or negative roots of C 2 . Note that both 
helds have opposite behavior for both roots while the behavior of arguments will be 
reversed if we take negative value of the root. The case |c 4 | < jcsl, gives imaginary form 
of the scale factor which was considered to be real. Hence we do not include it in our 
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analysis. 

Ic. (cj = 0 for all z = 2,5, 144ci — 3Cj = 0) 

a{t) = {Cit + , 

a(t) = const., 

(j)(t) = const., 


(83) 

(84) 

(85) 


which is a solution without dynamical helds. Note that for C 2 = 0, we obtain 
a oc which is the standard Friedman solution. Furthermore the constraint on ci, 
ensures that all invariants are also satished. Since Ji corresponds to time-translational 
symmetry with constant gauge, therefore Ji = Ci. The choice of this constant requires 
more understanding. In classical mechanics where a physical system is described by a 
Lagrangian that is the difference of a kinetic term and a potential term, the invariance 
of Lagrangian under time translation symmetry guarantee the conservation of energy 
(Hamiltonian). In this case we are not sure what this amounts to in our case as the 
given Lagrangian flSTl) does not carry the same interpretation as in classical physics. On 
the other hand the Hamiltonian is always assumed to be zero in gravity and general 
relativity lacks to provide a precise understanding of the concept of energy of the 
universe. Without loss of generality we can redehne Ci in terms of Ci, say Ci = ±Ci, 
where ci = 30^/144. Hence we get 

a{t) = {±cit + , ( 86 ) 


which gives us the same graph in either case (±), depending on the choice if t is assumed 
to be positive or negative. The positive values of C 2 , provides a translational shift to the 
graph on the left and on the right for negative values of 6 * 2 . It is pointed that although 
the form of potentials in above case is same as the form of potential in [1] (Section 4), 
where they obtained an exponential form of the scale factor. It does not contradict our 
analysis. Because we can obtain the same form by taking the value ci = 0, i.e. when the 
“Hamiltonian” (energy) is equated to zero. This case appears below for non-constant 
potential in a more general form. 

Case 2. (H(0, a) = C ^ 0) 

In this case the scaling symmetry evaporates and we get same four Noether symmetries 
Xi, where i = 1, 3,4, 5 and invariants. We hrst equate them to constants, = Cj, where 
z = 1, 3,4, 5. Now in this case the invariants J 4 and I 5 , implies that 


) *^4 . C5 

0=^, ^ = (87) 

O'* 

which reduces J 3 = C 3 , to a form that establishes a linear relationship between (j) and a, 
given by 040 " -|- c^cj) -|- C 3 = 0. The remaining invariant Ji, results into 


o 4-2 

3a a 


Ca^ 


= Cl, 


( 88 ) 


which can be simplihed into 


/c + Ca^y^^ 

V 3a4 ) 



a = 


C= Cl + 


2 


(89) 
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It can not be explicitly integrated. However, there arise a case in which if all gauge 
constants are zero c = 0 , then it results into following solution. 

2a. (cj = 0 for all i = 1, 3,4, 5) 


a{t) = exp - <^ 1 ) j > 

a{t) = const., 

0 (f) = const.. 


(90) 

(91) 

(92) 


where it is clear that potential {V = C) can not be negative and Ci is an integration 
constant which produces a shift in the graph on left if it is positive and right if it is 
negative. 

Case 3. (V (0, a) = (0 ± cr)", n 7 ^ 0) 

For this case, the Lagrangian fl54l) admits three Noether symmetries fl5^ . The action 
of the three dimensional group on the Lagrangian fl5T)) . in the form of continuous 
transformations of constant parameters, is given below 


Xi = dt —L(f - ei,a,0,a), 

X2 = d^Tda —t L(f,a, 0±e2,cT=Fe2), 

= tdt + ^da =F - (o-d^ + (t)da) —)• 0, d), 

3 n 


where, 0 = cosh f —- 
\ n 

/2e3 


= 1 = sinh (— ) a. 


\ n 


a = =F sinh 


\ n 


cf) + cosh 


n 


a. 


Therefore scaling symmetry (which is a consequence of homogeneity) arise with the 
requirement that both helds 0 and a, undergo a hyperbolic rotation such that the 
difference of their squares 0^ — cr^, is preserved. It is emphasized that the dynamical 
invariants in both cases are crucial to investigate quintom models. Therefore we hrst 
bring them in a more suitable form to apply further analysis. It is interesting to see 
that Is, in case 2, can be simplihed into a form 


2 ar 

I 3 = 2 a^a =F-(0d — ct 0) — f/i, 

n 


(93) 


using the dehnition of Ji. Therefore we start with the assumption ci 7 ^ 0, so we obtain 


2 a'^ 

Is = 2a^a =F-( 0 d — cr(j)) — cit = C 3 , 

n 


(94) 


where C 3 7 ^ 0. In Cosmology, the equations are generally described in Hubble parameter, 
H = a/a, in which case the above equation takes the form 


2 a-^ 


H =F — (0CT — CT0) 
n 


Cit = C3. 


(95) 


The other invariant I 2 is 

Jo = ( 0 ± d 


= C 2 . 


(96) 
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Here we arrive at following subcases. 

3a. (c 2 = 0) 

Note that C 2 = 0, implies that (p = =Fcr, which can be integrated to get 0 = + Ci, 

where Ci is a constant of integration. Therefore equation (l95|) reduces to 

H =F - cit = C 3 . (97) 

n J 



This is a constraint equation on the equations of motions (|2]), as provided by the Noether 
gauge symmetry and help to reduce the dynamics. Note that it is a first order ordinary 
differential equation in two unknown functions and we require Ci 7 ^ 0 , as otherwise 
(Cl = 0) the Lagrangian ([H]) becomes trivial. We solve equations of motions (E]), 
along with above constraint equation and obtain following cases where evolution of 
all dynamical quantities is also given. 

3a. (i) 1/ = (0 + a)"'. 

I- {n 7 ^ 1, Cl = 0, C 3 = 0) 


ait) 

a{t) 




± 




(n-2)' 


1/n 


nt + C3, 


II- (n 7 ^ 1, Cl = 0, C 3 7 ^ 0) 

a(() = 

. t 

a{t) = 


ncs e^ 


■2C|(3C"+2)i/2 


± 


C" 


n—2 


3a. (ii) V = {(f) — cr)"' 

I- (n 7 ^ 1, Cl = 0, C 3 = 0) 


1/2 

nt + C3, 


(98) 

(99) 


( 100 ) 

( 101 ) 


a{t) = •, 


a(t} = q: 


C! 


(n-2)' 


1/n 


nt + C3, 


II- (n 7 ^ 1, Cl = 0, C 3 7 ^ 0) 

a(() = C2e±vVT73) 

t 


ait) = T 


nc3 e' 


2C|(3C”+")V2 


T 


Cl 


- 2 \ 1/2 


nt + C3, 


( 102 ) 

(103) 


(104) 

(105) 


3b. (c 2 7 ^ 0) 

For this case we were unable to find any closed form of the solutions, however, invariants 
and provide us constraint equations to be satisfied by any quintom model. We 
have also verihed that any possible solution to these equations yield a viable quintom 
model. This is mainly done by solving above equations (hve of them) along with 
the condition that tc' = 0 , from equation ([7]), it turned out that the entire system 
is inconsistent. Therefore all such quintom models require that tc' 7 ^ 0. 













Noether symmetry in quintom cosmology 


14 


Case 4. (V (0, a) = (0 =F cr)"'/ (0 ± a)™, n 7 ^ 0 7 ^ m) 

For the last case we first write invariants fl641) . in a proper form. For example we can 
use Ji = Cl, relative to time translation to simplify I 2 , in (IMll . where we get 



(106) 


a n + m 


Again in this case we were unable to hnd closed form of the solutions. However, equation 
01061) . provide us a constraint equation to be satished by any viable quintom model. 
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